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5 s '. Abstract 

We study the orbital magnetism of 2D lattices with chaotic motion of elec- 
trons withing a primitive cell. Using the temperature diagrammatic technique 
we evaluate the averaged value and rms fluctuation of magnetic response in 
the diffusive regime withing the model of non-interacting electrons. The fluc- 
tuations of magnetic susceptibility turn out to be large and at low temperature 
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■ can be of the order of xl(&f0 3 ^ 2 ) where kp is the Fermi wavevector, I is the 
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mean free path, and xl is the Landau susceptibility. In the certain region 



of magnetic fields the paramagnetic contribution to the averaged response is 
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field independent and larger than the absolute value of Landau response. 



In the recent years, the orbital magnetism of mesoscopic systems has attracted much 

•rH . 

attention ||T]-§J. In the limit of the classically weak magnetic field, when the cyclotron radius 
is larger than the electron mean free path, the structure of energy levels of mesoscopic sample 
is very sensitive to the impurity configuration. As a result, the fluctuations of magnetic 
response turn out to be larger than the disorder-averaged value, i.e., the magnetic response 
of mesoscopic metallic sample has a random sign as a function of impurity configuration. 
This makes the investigation of sample-specific fluctuations to be of great importance for 
the understanding of the orbital magnetism of mesoscopic systems. 

In the work of Oh et al. the most attention was paid to the extreme quantum coherence 
limit, L <C Lh, Lt, L^. Here L is the linear size of sample, is the phase coherent length, 



Lh — \ ch/eH is the magnetic length, = JhD/T is the thermal length, D is diffusion 



coefficient and T is the temperature. Mesoscopic fluctuations in the limit Lt < L were 
calculated by Raveh and Shapiro The averaged magnetic response in different regimes 
was considered by Altshuler et al. ||. The authors of these works dealt with the case of 
random distribution of impurities all over the system. 

In this letter, we study the orbital magnetism in 2D chaotic lattices, i.e., in periodic 
complex systems with chaotic motion of electrons withing a primitive cell. As example of 
such systems, we consider 2D metallic sample (linear size L) consists of N parts (linear size 
a = L/\fN) with the same random configuration of impurities ||. We assume that a»l, 
so that the motion of electrons withing a primitive cell is diffusive and can be described by 
the diffusion constant D. The system under consideration is schematically shown in Fig. 1. 

We evaluate the mesoscopic fluctuations of magnetic response from the correlation func- 
tion of the thermodynamic potential Q. The latter is related to the correlation function of 
particles number M by the following expression: 



where 8J\fi = 8Af(p>i; Hi), p is the chemical potential, H is the magnetic field. The symbol 
( )i a t means averaging over the different configurations of impurities which have the long- 
ranged correlation of their positions. 

The correlation function of particles number is related to the correlation function of 
electron density p(r) : (8Ni8N 2 )i a t = Is Is dridr 2 (8pi(r 1 )8p 2 (r2))iat, where S = L 2 is the 
square of lattice. Scattering contribution to the electron density at some point r t is due 
to the sum over the all possible scattering paths returning to ri. The correlation function 
(8pi(ri)8p2(r 2 ))iat is determined only by the contribution of the pairs of phase-coherent 
electron paths (see, e.g. Qj). In the case of lattice, an electron traversing the closed path 
starting at the point ri and passing near the point r 2 + N acquires the same phase as an 
electron traversing the similar, but shifted in space by N, path, as it is shown in Fig. 1. 
Here N is 2D vector with components {ai,aj}, where integers i,j = 0...\/~N — 1 numerate 
the primitive cells of lattice. Thus, calculating the sum over N we obtain 
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(*«> w = iV(MW 2 )o , (2) 

where ( )g means averaging over the paths shifted by N = 0. 

The correlation function (5J\fi8J\f 2 )fi can be conveniently calculated by the temperature 
diagrammatic technique. We use the diffusion approximation which is valid when T > A, 
where A = (a 2 i>)~~ 1 is the mean spacing between energy bands in chaotic lattice, v being the 
mean density of states HH- The largest contributions to (AA/i&A/^o are the double-diffuson 
and double-Cooperon diagrams. We consider magnetic fields applied perpendicularly to 2D 
system and satisfying the condition I <C Lh- We also assume that Lt < L. For this case 
the correlation function is 

ST 1/r / 1 \~ 2 

{8N X 8N 2 )$ = —— J2uReY,Y. UJ 4 UJ + huJ s( n + 9) + ^12 • (3) 

n ~ nIJ u»0 n>0s=± V I J 

Here u = 2nTm, m is an integer number, /ii 2 = \i\ — fi2, oj± = 2De(H\ ± H 2 )/(ch), 
c is the speed of light, e is the electron charge, n is the Landau level number. The upper 
cutoff in (3), 1/r, is due to the validity of the diffusion regime. After substitution of (2) and 
(3) in (1) and some transformations we arrive at the following expression for the correlation 
function of thermodynamic potential in the form of t-dependent integral: 

(6Qm6Q(H 2 )) lat = J q t2sinh2(7rTt) [ sinh{uJ+t/2) + sinHu _ t/2) ) ■ ( 4 ) 

To obtain the correlation function of the magnetic moment, one has to successively 
differentiate (4) on Hi and H 2 . We would like to note that the integral in (4) converges only 
after differentiation on Hi and H 2 . The analytical form of the correlation function can be 
found in two limiting cases: 

1) In the "high temperature" limit, huo + T, we obtain 
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(5M(Hi)5M(H 2 )) lat = - — (k F l) 2 X 2 L HiH 2 . (5) 

Here xl = —e 2 S/(12Trm*c 2 ) is 2D Landau susceptibility, m* is an effective electron mass. 
The rms fluctuation of magnetic susceptibility can be directly obtained from (5): 
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(Sx 2 )ll?= 2 ^(k F l)\XL\. (6) 

If impurity positions are random withing entire system (no long-range correlation, i.e., N = 
1), a should be substituted on L. In this case, Eq. (6) coincides (up to a numerical factor) 
with the result of Raveh and Shapiro [[J. Here we would like to draw attention to the 
following peculiarity of chaotic lattice. We remind that Eq. (6) was obtained for the case 
Lt/L < 1. At the same time, the quantity L F /a = NL F /L, appearing in the Eq. (6), 
can be large. Indeed, the validity of the perturbative approach is T > A and in the low 
temperature limit T ~ A we have f\L T /a ~ (HD/Aa 2 ) 1/2 ~ (hDu) 1/2 ~ (k F l) 1/2 . Therefore, 
the rms fluctuation of magnetic susceptibility in the case of lattice can be of the order of 

(M) 3/2 |xl|. 

2) In the "high field" limit, Tiu + ^> T, we have 

(5M(F 1 )5M(iJ 2 ))^ ot = l^l^(M)\i(|^i + H 2 \ - \H X - H 2 \) , (7) 

where ((x) is Riemann Zeta- Function. Such form of the correlation function indicates that 
the magnetic moment as well as the amplitude of its oscillations are growing functions of mag- 
netic field H. It is also interesting to note that the correlation function (5 M (Hi)5 M (H 2 )) i a t 
depends only on min{Hi, H 2 }. 

Now let us calculate the canonically averaged magnetic response (M)i at . It can be 
expressed as the sum of the diamagnetic Landau response and a paramagnetic contribution: 
(M)i at = X lH + (Mp)fat, where 

(M p ) lat = -^^(5Ar 2 ) lat . (8) 

For complete discussion about the difference between canonical and grand-canonical situa- 
tion see and references therein. Note that in the case of lattice, the number of electrons 



1 We can consider the temperatures T ~ A when it does not contradict to the condition Lt < L 
(T > A(k F l)/N), i.e., for the lattices with N > (k F l) > 1. 
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in a primitive cell is fixed. Calculating the correlation function (5N 2 )i a t (Eqs. (2) and (3) 
at H\ = /i 2 = and H 1 = H 2 = H) we finally obtain 



±^(k F l)\ XL \H , if ^«T 

W to = <; (9) 

£MeM hu H »T. 

Here tu// = ADeH/ (ch) is the Cooperon cyclotron frequency. In the weak magnetic fields, 
hujH <C 7" (-^ij S> -^t) the averaged magnetic response is paramagnetic and grows linearly 
with H. In stronger magnetic fields, Tiujh ^ 7" {Lh "C -^t) ; the paramagnetic contribution 
(Mp)i at is H - independent. As one can see from (9), the paramagnetic contribution (M p )i at 
is larger than the absolute value of diamagnetic Landau response \xlH\ until magnetic 
fields hujn ~ A(kpl) (Lh ~ a). We note that in the case of small metallic particles the 
paramagnetic contribution (M p ) becomes if-independent only in the fields, where it is small 
compared to the Landau response j3|. 

In conclusion, we discuss the approximation that we used in derivation of (3). In the 
calculation of diffusion propagators we neglected a contribution from another type of phase- 
coherent electron paths (fragment of such trajectories shown in Fig. 2). The contribution 
from the paths shifted by 2N is of the order of weak localization correction to the diffusion 
constant 

5D(2N) 1 TJ ^ T . 

where W(R,uj) is the probability to diffuse on the distance R during a time u~ x . Be- 
cause the number of such trajectories is of the order of (Lt/o) 2 , the total contribution 
^ N SD{2N)/D ~ (L T /afl/(k F l) ~ A/T. For T > A we can neglect it. We note that this 
correction depends on the magnetic field. This dependence contains not only a monotonic 
part, but also a part which oscillates with the period $ c / $o, where $ c = Ho 2 is the magnetic 
flux through a unit primitive cell, <3>o — hc/e. 
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